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I.  INTRODUCTION 


For  almost  50  years  there  has  been  great  experimental*  and  theo¬ 


retical' 


interest  in  the  gas-phase  recombination  of  hydrogen,  the  process 


H  +  H  +  M  -*■  H2(v,J)  +  M  (1) 

Reaction  (1),  where  M  =  H,  represents  one  of  the  simplest  chemical  reactions 
and  is  of  fundamental  interest  to  chemical  dynamics  and  astrophysics.  Indeed, 
this  reaction  tests  the  outer  limits  of  the  (H)j  potential  energy  surface. 

This  reaction  is  of  particular  interest  because  of  the  recent  calculations  and 
experiments  on  the  similar  H  +  D2  system.  The  comparison  of  theory  and 

experiment  in  these  simple  three-electron  systems  tests  our  ability  to  under¬ 
stand  chemical  reactions. 

Recombination  of  hydrogen  is  of  practical  importance  in  both  plasaia 
physics  and  rocket  propulsion,  especially  for  hybrid  and  electrothermal  pro¬ 
pulsion  systems .^*2^  This  reaction  becomes  important  in  chemical  rockets  if 
hydrogen  is  a  major  atomic  component  of  the  fuel.  In  these  systems,  the  total 
recombination  rate  coefficient  and  the  distribution  of  energy  into  the  product 
vibrational /rotational  states  are  important.  For  efficient  propulsion,  the 
translational  energy  of  the  exhaust  must  be  maximized.  Frozen  flow  efficiency 
losses  occur  when  the  molecules  do  not  sustain  enough  collisions  in  the  nozzle 
to  equilibrate  the  translational  and  vibrational/rotational  degrees  of  free¬ 
dom.  The  energy  that  is  trapped  in  the  vibrational/rotational  modes  of  the 
system  is  not  available  for  thrust. 

Experimental  reaction  rate  data  were  critically  reviewed  by  Baulch  et  al. 
in  19722^  and  by  Cohen  and  Westberg  in  1982. For  many  third  bodies,  includ¬ 
ing  M  *  H2,  He,  and  Ar,  accurate  measurements^*^  and  accurate  calculations* 2 
are  available  and  are  in  reasonable  agreement.  For  M  ■  H,  there  is  a  lack  of 

reliable  experimental  data.  Shock  tube  measurements  at  high  temperatures  have 
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been  performed  but  have  not  been  very  reproducible.  *  At  lower  tempera¬ 
tures,  only  an  upper  limit  to  the  rate  coefficient  at  300  K  is  available, 
which  was  obtained  by  Bennett  and  Blackmote.  This  upper  limit  is  net  in 
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agreement  with  the  theoretical  result,12’17  and  the  experimental  method  has 
been  questioned  by  Baulch  et  al.2^ 

The  nascent  vibrational/rotational  product  state  distribution  of  recom¬ 
bined  H2  has  never  been  measured.  In  this  report,  we  present  a  calculation  of 
this  distribution  using  resonance  complex  theory.  The  theory  and  computa¬ 
tional  details  are  described  in  Section  II.  The  total  recombination  rate 
coefficient  is  presented  and  compared  to  previous  calculations  in  Section 
III.  The  nascent  vibrational/rotational  distribution  is  discussed,  and  rate 
coefficients  into  individual  vibrational /rotational  states  as  functions  of 
temperature  are  presented. 


II.  COMPUTATIONAL  DETAILS 

A  three-body  recombination  process  can  take  place  primarily  via  two 
separate  channels:  the  "energy-transfer”  mechanism 

H  +  H  ^  H2* 


H2  +  M  -*■  H2  +  M 


and  the  "exchange"  or  "chaperone"  mechanism 

H  +  M  ^  HM*  (3a) 

HM*  +  H  H2  +  M  (3b) 

where  H2  and  HM  are  intermediate  (unbound)  complexes  that  are  in  equil¬ 
ibrium.  In  the  present  case,  where  M  =  H,  the  two  processes  can  be  distin¬ 
guished  in  a  classical  calculation;  quantum  mechanically  these  two  channels 
are  indistinguishable. 

The  resonance  complex  theory  of  Roberts,  Bernstein,  and  Curtiss**  is  used 

to  study  the  recombination  process.  Because  this  theory  has  been  discussed  in 

detail  elsewhere,  only  a  brief  description  follows.  In  this  theory,  the 
complex  [H2*  in  Eq.  (2)]  is  Identified  with  orbiting  resonance  states  of  the 
H2  system.  These  states  are  bound  classically  but  have  a  finite  lifetime 
quantum  mechanically,  because  of  tunneling.  We  assume  that  some  of  these 
states  reach  a  steady-state  population  distribution  determined  by  their 
equilibrium  constant.  Hence  the  reaction  rate  can  be  written  as  a  sum  over 
Individual  rate  constants 

kr(v,J,T)  k*(v,J,T)-  (4) 

i 

The  individual  rate  coefficient  for  deactivation  of  an  orbiting  resonance  i 


can  be  expressed  as 


4<v,J,T)  -  ;rK‘  J,T) 


B8 


ten 
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(5) 
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where 


=  (8kT/np) 


K*  is  the  equilibrium  constant  for  the  formation  of  the  ic^  resonance  complex 
state,  u  is  the  reduced  mass  for  the  (l^M)  system,  and 

00 

oi(v,J,T)  =  (kT)~2  /  E0i(v,J,E)exp(-E/kT)dE  (7) 

o 

is  the  thermal  average  of  the  energy  dependent,  state-selective,  cross  sec¬ 
tion. 

To  calculate  vibrational/rotational  rate  coefficients,  the  resonance 

states  contributing  to  the  hydrogen  (H  +  H  +  H)  recombination  process  must  be 
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determined.  There  are  approximately  50  orbiting  resonance  states  of  H 2 .  » 

The  procedure  for  assessing  the  contribution  of  each  resonance  state  was 

discussed  in  the  first  calculation  using  resonance  complex  theory.11  Most  of 

these  states  can  be  eliminated  by  simple  energy  considerations  because  the 

equilibrium  constant  contains  a  factor  of  exp(-E/kT).  The  majority  of  the 

remaining  states  (about  15)  have  lifetimes  mich  too  long  (i.e.  ,  too  narrow  in 

resonance  width)  to  contribute  to  the  recombination.  The  result  is  6  possible 

contributing  orbiting  resonance  states.  The  properties  of  these  states  are 

summarized  in  Table  1.  The  literature  includes  discussions  about  which  of 

these  6  resonance  states  contributes  to  recombination  under  various  experi- 
6  8 

mental  conditions.  ’  We  decided  to  calculate  contributions  from  all  6 
states.  If  it  is  determined  that,  under  certain  experimental  conditions,  some 
of  the  states  do  not  contribute  or  are  not  in  complete  thermal  equilibrium, 
they  could  then  easily  be  excluded  from  the  final  calculation  of  the  rate 
coefficient . 


raw** 


Table  1.  Properties  of  Six  Resonance  States  in  Recombination  Calculation 


(cm" ^ ) 


<R>(au) 


T^sec) 


14 

5 

45.7 

6.51 

2.7  x  1 

14 

4 

3.6 

6.05 

4.0  x  1 

13 

8 

89.8 

5.31 

2.0  x  1 

12 

1 1 

216.0 

4.86 

1.8  x  1 

12 

12 

387.0 

5.30 

7.0  x  1 

1 1 

13 

199-5 

4.27 

2.0  x  1 

Two  of  the  resonance  states  require  an  additional  approximation.  The  v  = 
1A,  J  =  5  and  the  v  =  12,  J  =  12  states  are  classically  unbound.  In  order  to 
carry  out  classical  trajectories  with  these  initial  states,  it  is  therefore 
necessary  to  select  them  to  be  slightly  classically  bound,  that  is,  with  a 

29 

slightly  lower  internal  energy.  The  vibrational  quantum  number  defined  by 


v 


-1/2  +  (1/uTO  J 


f  2  1/2 

{2m[Eint  -  V (R)  -  (JrJr/2mR  )]}  dR 


(8) 


where  are  the  diatomic  turning  points,  E^nt  is  the  internal  energy,  V(R)  is 
the  diatomic  potential  energy,  m  is  the  diatomic  reduced  mass,  and  Jr  is 
related  to  the  rotational  quantum  number  by 


=  j(j  +  l)h2 


(9) 


then  takes  on  noninteger  values.  The  "vibrational  quantum  number"  for  the 
v  =  1A  resonance  state  is  v  =  13.85  and  for  the  v  =  12  state  is  v  =  11.72.  To 
assess  the  effects  of  this  change,  calculations  are  run  for  not  only  these  two 
states  but  also  for  other  nearby  (E^nt  slightly  different)  states.  No  effect 
is  found  either  on  the  total  cross  section  or  on  the  final  vibrational/ 
rotational  state  distribution. 

Once  the  resonance  state  is  selected,  quasiclassical  trajectories  are 
calculated  for  the  process 

H2*(v',j')  +  H  -*■  H+H+H  (10a) 

♦  H2(v,J)  +  H  (10b) 


to  converge  o^(v,J,E).  A  similar  calculation  with  M  *  H,  He,  Ar,  and  H2  was 
performed  by  Whitlock  et  al.*2  for  the  total  recombination  rate  coefficient. 
The  present  calculation  follows  Che  same  basic  procedure.  However,  several 
changes  were  needed  to  determine  the  vibrational/rotational  distribution  of 
products.  Classical  trajectories  on  the  (H)^  potential  energy  surface  used  by 
Whitlock  et  al . were  calculated  using  CLASTR.^®  We  discovered  that  the 
surface  supports  many  bound  vibrational/rotational  states  that  are  unbound  for 
the  true  H2  system.  A  highly  accurate  surface  for  H-j  is  available,  the 

12 


v' V -V vV 'v' V V v V v -v 


SLTHJi»  surface.  Therefore,  we  modified  CLASTR  to  accept  this  surface  and 
to  determine  the  final  vibrational  and  rotational  state  quantum  number  for  the 
SLTH  surface.  We  constructed  the  correct  initial  conditions  for  the  tra¬ 
jectories  based  on  the  orbiting  resonance  state  and  calculated  the  final 

vibrational  and  rotational  quantum  numbers.  The  procedure  is  similar  to  the 

2  9 

one  described  by  Truhlar  and  Muckerman  and  will  not  be  reviewed.  Because  of 
the  highly  excited  states  encountered  for  both  the  initial  and  final  condi¬ 
tions,  we  made  one  modification  to  their  described  procedure.  The  Newton- 
Raphson  iteration  technique,  to  locate  the  turning  points  of  the  effective 
potential,  was  replaced  with  the  method  of  Pade  approximates. 

After  the  individual  vibrational /rotational  cross  section  has  been  cal¬ 
culated  as  a  function  of  E,  it  is  then  fitted  to  the  form 


o ^ (E  , v ,  J )  =  b(E/k)-1/2exp[  -a(E/k)-1] 


as  suggested  by  Whitlock  et  al.  After  thermal  averaging  to  determine 
o^Cv.J.T),  the  cross  sections  are  combined  to  produce  the  rate  coefficient 
for  the  desired  vibrational/rotational  state. 

For  each  of  the  six  resonance  states,  five  energies  are  calculated  (E/k  = 
50,  100,  300,  1000,  and  2000  K).  This  is  found  to  be  adequate  to  calculate 
the  thermally  averaged  cross  section.  For  each  energy  and  initial  resonance 
state,  approximately  1000  trajectories  were  run,  producing  a  Monte  Carlo  error 
of  less  than  5 %  for  the  total  cross  section,  of  approximately  10%  for  the 
summed  vibrational  cross  sections,  and  of  approximately  20%  for  the  individual 
vibrational/rotatational  cross  sections. 


III.  RESULTS 


The  individual  final  state  vibrational/rotational  probabilities  for  one 
initial  resonance  state,  v  =  13,  J  ■  8,  are  shown  for  E/k  *  50  K  in  Fig.  1  and 
for  E/k  *  1000  K  in  Fig.  2.  The  qualitative  features  are  the  same  for  all  the 
resonance  states.  The  dissociation  level  for  the  hydrogen  molecule  is  shown 
by  the  heavy  line  in  the  figures.  The  figures  indicate  that  the  distribution 
is  peaked  near  this  line,  that  is,  in  the  highest  bound  states  of  the  H2 
manifold.  The  distribution  broadens  and  extends  to  lower  vibrational  states 
as  the  collision  energy  increases.  The  total  reaction  probability  decreases 
with  increasing  energy.  Figure  3  exhibits  the  individual  vibrational  cross 
sections  for  the  same  intermediate  state  as  functions  of  collision  energy. 
Again,  the  qualitative  features  are  the  same  for  all  resonance  states.  In 
this  figure,  the  widening  of  the  vibrational  distribution  and  the  f alloff  of 
the  total  cross  section  with  increasing  collision  energy  can  be  seen 
clearly.  (Space  considerations  preclude  the  display  of  all  the  individual 
vibrational /rotational  cross  sections  as  functions  of  energy.  The  complete 
results  are  available  on  request.)  The  nascent  vibrational/rotational 
distribution  is  the  correctly  weighted  sum  over  the  individual  orbiting 
resonance  states. 

In  Fig.  4,  the  results  for  the  total  rate  coefficient  as  a  function  of 
temperature  are  plotted  and  compared  with  the  calculation  of  Whitlock 
et  al.*^  The  two  data  sets  reveal  the  6ame  qualitative  behavior  with  tempera¬ 
ture.  However,  the  present  results  are  approximately  a  factor  of  2  smaller 
than  the  previous  results.  This  difference  occurs  because  the  potential 
energy  surface  used  by  Whitlock  et  al.^  binds  too  many  vibrational /rotational 
states,  as  discussed  previously.  In  Fig.  5,  the  rate  coefficients  into  indi¬ 
vidual  vibrational  states  are  plotted  and  compared  with  the  total  rate  coeffi¬ 
cient  and  with  the  sum  over  the  plotted  states.  This  figure  indicates  that 
v  *  14,  v  *  13,  v  ■  12,  and  v  ■  11  make  up  the  largest  contribution  to  the 
rate  coefficient.  Figure  5  also  exhibits  the  broadening  of  the  distribution 
with  energy,  with  v  *  14  dominating  the  rate  constant  at  lower  temperatures 
but  becoming  overshadowed  by  v  =  13  at  higher  temperatures.  Figure  6  displays 
the  rotational  rate  coefficients  for  v  ■  14.  The  higher  rotational  levels 
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Fig.  A.  Total  Recombination  Rate  Constant  as  a  Function  of  Temperature 


TOTAL  RATE  CONSTANT 

SUM  OF  RATE  CONSTANTS  INTO  FOUR  HIGHEST 
LYING  VIBRATIONAL  STATES 


dominate;  J  *  3  is  the  highest  bound  state  for  v  *  14.  Figure  7  exhibits  the 
rotational  rate  coefficients  for  v  =  13.  Again,  the  higher  rotational  states 
dominate,  with  the  highest  four  states  constituting  approximately  80%  of  the 
total  rate  coefficients.  (J  *  7  is  the  highest  bound  state  for  v  *  13.) 

These  results  indicate  that  the  nascent  vibrational/rotational  distribu¬ 
tion  peaks  in  the  highest  bound  states  of  the  hydrogen  molecule.  The  average 
vibrational  level  is  high.  Because  such  high  vibrational  states  can  only  sup¬ 
port  a  small  number  of  bound  rotational  states,  the  average  rotational  level 
is  small. 


IV.  CONCLUSIONS 


We  have  used  resonance  complex  theory  to  investigate  the  nascent  vibra¬ 
tional/rotational  state  distribution  produced  by  hydrogen  atom  recombina¬ 
tion.  We  have  shown  that  this  distribution  peaks  in  the  highest  bound  states 
of  the  hydrogen  molecule.  Because  of  the  highly  excited  nature  of  the  pro¬ 
ducts,  the  nascent  distribution  is  greatly  changed  by  subsequent  collisions. 

To  assess  the  effect  of  a  collision  on  the  distribution,  calculations  must  be 
performed  to  determine  the  rate  constant  for  dissociation  and  relaxation  from 
these  states.  This  is  particularly  important  when  inert  third  bodies  such  as 
helium  are  considered.  The  measured  hydrogen  ortho/para  ratio  in  such  systems 
may  not  be  the  nascent  distribution  because  of  the  effect  of  vibrationally 
enhanced  exchange  reactions. 

This  state  distribution  was  calculated  using  the  approximations  of  quasi- 
classical  theory.  Though  a  correct  quant umnmechani cal  state  was  used  for 
initial  conditions,  the  subsequent  collision  was  treated  using  classical 
mechanics.  As  has  been  previously  discussed,  classical  results  are  more 

O  Q 

accurate  the  greater  the  degree  of  averaging  in  the  calculation.  Therefore 
we  would  expect  the  total  rate  coefficient  to  have  the  greatest  accuracy,  the 
vibrational  rate  constants  to  be  less  accurate,  and  the  individual 
vibrational /rotational  rate  constants  to  be  the  least  accurate.  The  qualita¬ 
tive  features  of  the  distribution  should  be  correct,  that  is,  the  distribution 
is  almost  certainly  peaked  in  the  highest  states  of  the  bound  hydrogen  mole¬ 
cule. 

Future  calculations  will  address  the  effects  of  a  collision  on  the  nas¬ 
cent  distributions  and  will  investigate  the  nascent  vibrational /rotational 
distribution  produced  with  an  inert  third  body  such  as  helium.  We  hope  that 
the  results  will  be  more  readily  subject  to  experimental  test. 
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The  Aerospace  Corporation  functions  as  an  "architect -engineer"  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corpora! iou ' s  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff’s  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  svstems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out -of -f ield-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  faul t -tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratoty:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  microwave/mi 1 1 imeter  wave  measurements, 
diagnostics  and  radiometry,  microwave/mil 1 imeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas ,  rf  systems,  electromagnet ic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  **11  ir.  spare  and  enemy -■* nd need 
environments . 

Space  Sciences  Laboratory:  Magnetospher ic ,  auroral  and  cosmic  ray 
physics,  wave-particle  Interactions,  magnetospheric  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
i  ns  t  r  urne  n  t  a  t  i  o  n . 


